Interplay between charge-order, ferroelectricity and ferroelasticity: 
tungsten bronze structures as a playground for multiferroicity 
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Large electron-electron Coulomb-interactions in correlated systems can lead to a periodic arrange- 
ment of localized electrons, the so called charge- order^\ The latter is here proposed as a driving 
force behind ferroelectricity in iron fluoride Ko.eFeJj'eFefj'iFs. By means of density functional the- 
ory, we propose different non-centrosymmetric d^/d^ charge-ordering patterns, each giving rise to 
polarization along different crystallographic axes and with different magnitudes. Accordingly, we 
introduce the concept of "ferroelectric anistropy" (pecuUar to improper ferroelectrics with polar- 
ization induced by electronic degrees of freedom), denoting the small energy difference between 
competing charge- ordered states that might be stabilized upon electrical field-cooling. Moreover, 
we suggest a novel type of charge- order-induced ferroelasticity: first-principles simulations predict 
a monocUnic distortion to be driven by a specific charge-ordering pattern, which, in turn, unam- 
biguously determines the direction of ferroelectric polarization. Ko.eFeo^eFeo^iFs therefore emerges 
as a prototypical compound, in which the intimately coupled electronic and structural degrees of 
freedom result in a manifest and peculiar multiferroicity. 

PACS numbers: Valid PACS appear here 



Materials which combine magnetism and ferroelectric- 
ity, belonging to the intriguing class of multiferroics, can 
be classified into two categories: [Ij "Structural magnetic 
ferroelectrics" (SMF), where the primary order param- 
eter in the ferroelectric (FE) phase transition is related 
to a structural instability (which can be polar or non- 
polar, see for example the prototypical BiFeOs or the 
"geometric FE" YMnOs), and "electronic magnetic fer- 
roelectrics" (EMFs), where the primary order parameter 
is related to electronic degrees of freedom, such as spin, 
charge or orbital order. [2j, |3l, 0, [5[ 

Whereas plenty of studies have been recently per- 
formed in the field of spin-driven ferroelectricity, ferro- 
electricity induced by charge-order (CO) still constitutes 
a largely unexplored territory. Even in the two paradig- 
matic cases in which the pattern of Fe^+ and Fe^+ ions 
was suggested to break space inversion symmetry, e.g., 
LuFe2 04[i and Fe3O4 0,ll, the actual occurrence of fer- 
roelectricity seems controversial. 

Within this context, non-oxides - and fluorides in par- 
ticular - are interesting candidates as potentially novel 
improper multiferroics. Indeed, Ko.6FeQ^eFeo^4F3 which 
crystallizes in a non-centrosymmetric tetragonal tung- 
sten bronze (TTB) structure, has been reported to show 
CO, although the exact CO pattern is still debated. Q 
The family of TTB compounds, obtained by substi- 
tuting either tungsten with transition metals or oxy- 
gen with other anions, exhibit a lot of functionali- 
ties, such as ferroelectricity, ferroelasticity, pyro-/piezo- 
electric properties [lo|- Despite such a remarkable tech- 
nological appeal, the electronic structure of TTB com- 
pounds has not been deeply investigated from the the- 
ory point of view, mainly because of the complex crys- 
tal structure. Guided by the high potential of TTB 
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FIG. 1: C-centered primitive unit cell of the optimized crystal 
structure of Ba2NaNb50i5 in (a) ab plane and in (b) ac plane 
(only three Nb-Oe octahedrons are shown). Upon ionic relax- 
ation, all Nb (O) ions are displaced toward +z (-z) direction, 
so as to induce a large Pz. 
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multiferroics field, we here focus on 
Ko.6Feo^gFeo^4F3 showing ferroelasticity and ferroelec- 
tricity driven by non-centroymmetric Fe-d^/d^ CO pat- 
terns. To outline the novel physics in Ko.6Feo^6Feo^4F3 
we will first discuss the prototypical TTB compound, 
Ba2NaNb5 0i5, where ferroelectricity is driven by conven- 
tional Nb^+ off-centering and where no coupling between 
ferroelectricity and ferroelasticity is observed. 



Methodology and structural details - DFT simulatio] 
were performed using the VASP code and the PA^ 
pseudopotentials [l2f within the GGA+[/ formahsmp~1 
(U =5 eV and J=0 eV for Fe d-states). The cut-off ener^ 
for the plane-wave expansion of the wave-functions w; 
set to 400 eV and a k-point sheh of (2, 2, 4) was used f 
the Brihouin zone integration. The internal atomic coo 
dinates were fully optimized until the atomic forces we 
less than 0.01 eV/A while the lattice parameters we 
taken from experiments. The FE polarization P was Cc 
culated using the Berry phase method [14], by comparii 
the FE and paraelectric (PE) state; the latter is co: 
structed by imposing the z reflection in the atom 
structure. For simplicity, Ko.6FeQ^gFeQ^4F3 was treated j 
a ferromagnet; further complexity of the experimental 
suggested ferri- or weak ferro-magnetic spin configuratic 
is not expected to affect CO (and related ferroelectricit}' ^ , 
as in LuFe204.[i5[ 

d^-ness at Ba2 Na Nb^ 5 - At room temperature, 
Ba2NaNb50i5(BNN)[3 crystallizes in the polar or- 
thorhombic Cmm2 structure with a = 17.626, h = 
17.592, c = 3.995 A[17]. At low temperature, BNN 
shows quasi-commensurate and incommensurate phases 
with ferroelastic transition, leading to a crystal cell in the 
ab plane larger than the above mentioned unit celLp^sj 
Since ferroelectricity was reported not to be coupled with 
ferroelasticity|19|] and since our main focus is on FE prop- 
erties, we optimized BNN in the Cmm2 structure. As 
shown in FigUl the optimized NbOe octahedrons are sig- 
nificantly distorted from the tetragonal symmetry. The 
polar distortion along the z axis is driven by the off- 
center shift of Nb atoms, due to a strong hybridization 
between Nb empty Ad and O 2p states, with an average 
Nb^+ ionic displacement (with respect to the side O ions) 
of 0.16 A. The FE behaviour can be clearly interpreted 
on the basis of the "(i^-ness" criterion^^- The latter 
doesn't only cause electric polarization, but also results 
in a wide energy gap: 2.7eV in the FE phase, with respect 
to 2.3eV in the PE state. The calculated polarization is 
34.3/iC/cm^ along the c direction; to our best knowledge, 
this is the first theoretical estimate reported in the lit- 
erature and is in good agreement with the experimental 
value of 40/iC/cm2[2l|. 

Charge order at Kq^qF^^ qF^q^F^, - In iron fluoride, 
the complexity of the different phase transitions is still 
a matter of debate. Earlier experiments showed cou- 
pled ferroelectric/ferroelastic transitions to occur at 490 
K,[22[ along with CO. [2^ More recently, several transi- 
tions were reported: [P] a first structural transition around 
570K from tetragonal to orthorhombic, a second transi- 
tion at 490K where Fe^+/Fe^+ CO occurs, and a third 
transition around 290K to monoclinic structure, cou- 
pled with ferroelasticity. Experimentally, the orthorhom- 
bic Pha2 structure was reported, with a = 12.751, 
b = 12.660 and c = 7.975A[9|, the unit ceh includ- 
ing two TTB layers along c. From the symmetry point 
of view, the Pba2 group has four symmetry operations: 
{E, C2z, + (1^0), Gy + (1^0)} so that FE polarization 




FIG. 2: Down-spin charge of Fe t2g state, 1 eV below Ey^ for 
the Pba2 symmetry at (a) z ~ 1/4 and (b) 2; ~ 3/4 planes. 
Circled ions are mixed- valence Fe^'^^ ions. The symmetry 
operations which relate the four Fe^'^^ are also shown, (c) 
Full CO patterns starting from partial CO shown in (b) can 
be obtained by turning four Fe^'^^ sites into two Fe^^ and two 
sites at z ~ 3/4 plane. The positions of (blue) circles 
correspond to the one in (b). 



is in principle allowed to be induced only along the z 
axis. Given the d^ / d^ electronic configuration of Fe ions, 
Pz by "d^-like" hybridization is of course not expected 
here, so that alternative mechanisms, such as CO, should 
be invoked to explain ferroelectricity. As discussed in 
ref.[9], the Fe ions show CO in two FeO layers. How- 
ever, the CO is not "full" in the experimental Pba2 struc- 
ture: although 12 Fe ions are supposed to be 'Fe^^{d^) 
and 8 Fe ions to be Ye^^{d^) due to the stoichiometry 
(Ki2Fe^2Fe8^^F6o/cell), the Pba2 crystal structure gives 
10 Fe2+, 6 Fe^+ and 4 Fe2-^+ ions (cfr Fig|2l^a),(b)). As 
expected, the mixed-valence Fe ions lead to a metallic 
behaviour (not shown), kept even after the structural op- 
timization. In order to obtain an insulating state (obvi- 
ously needed for ferroelectricity to develop), the symme- 
try must be reduced so as to get a fully charge-ordered 
Fe sublattice. The existence of a larger supercell which 
may stabilize a complete CO was previously suggested; [9| 
however, this would require an extremely high and unaf- 
fordable computational load. Therefore, we limited our 
study to the Pba2 unit cell, artificially breaking symme- 
tries so as to obtain a "full" CO. Here we assumed three 
CO patterns, CO-I, CO-II and CO-HI, each keeping one 
of the three symmetry operations and breaking the other 
two among cr^; + (^0) and + (^^0), respectively 
(cfr Fig.2(c)). Whereas the induced polarization P^ is 
allowed by the prototype Pba2 crystal symmetry, Py and 
Px are additionally allowed in COII and COIII, respec- 



tively. 



TABLE I: Total energy difference (meV/Fe), energy gap (eV) 
and induced FE polarization calculated by Berry phase P^^^^^ 
and by dipoles P^^p^^^ (/xC/cm^) at partial CO pattern in 
experimental crystal structure, optimized structure keeping 
experimental symmetry, full CO patterns; CO-I, CO-II, CO- 
III (with given symmetry). 
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After ionic optimization at each CO-pattern, the 
charge separation between Fe^+ and Fe^+, calculated by 
integrating the charge density in lA atomic radius, is 
equal to 0.365 e~; this is rather large compared to the 
value of Fe304 (^0.2 e~) [S]. Indeed, this is consistent 
with the expected weaker Fe-F hybridization compared 
to Fe-0 and put forward fluorides as better candidates 
for CO - compared to oxides - where a larger charge- 
disproportionation can be achieved. 

CO-induced ferroelectricity - We found that, among 
the assumed CO patterns as well as the experimen- 
tal/optimized state, con is the most stable state with 
the largest energy gap (cfr Tab. [I]). Besides, COII 
shows a sizeable polarization p^erry _ _q 59 ^C/cm^. 
This can be roughly understood as induced by local elec- 
tric "point-charge" dipoles which connect Fe^+ and Fe^+ 
ions, as previously suggested for Fe304[8]. The electric 
dipoles are calculated considering only four Fe ions, i.e. 
those originally located on the mixed valence sites. In 
fact, the CO pattern of the other Fe ions are identical 
to the high-symmetry Pba2 CO pattern, so that their 
net contribution to P would cancel out. As shown in 
FigHb) and (c), P^^^p°^^ is calculated as -0.57 /iC/cm^ in 
good agreement with p^^rry. ^.j^^g shows that polariza- 
tion has a purely electronic origin and that ionic dis- 
placements play a minor role. Note that the sign of 
Py can be switched by exchanging Fe^+ and Fe^+ ions; 
under an applied electric field Ey^ charge shifts among 
these four Fe sites is indeed expected to occur. In an 
analogous way, we calculated Px at CO-III pattern as 
pdipoie ^ pBerry ^ — 5^C/cm^, a valuc much larger than 
P at CO-II. In this case, too, the agreement between 
pdipoie pBerry -g remarkable. Finally, we note that, 
as expected, the CO-induced polarization, though size- 
able, is much smaller than in Ba2NaNb5 0i5 where the 
ionic degrees of freedom are the source of ferroelectricity. 

Charge and Orbital Order - The CO arrangement 
also causes a change in the Fe-t2oOrbital ordering (00). 
When comparing Fig|2] and Fig|3l the occupied t2g or- 
bitals (both at z=l/A and z=3/4 planes) modify their 
shape/direction to avoid as much as possible the nearest 




FIG. 3: (a) Down-spin charge of Fe t2g states, 1 eV below 
Pf in the optimized structure with the ground-state COII-ctx 
pattern. Circles show the Fe sites where Fe^"^^ mixed- valence 
ions are originally located in the experimental structure, (b) 
Electric dipoles arising among circled Fe sites, due to the dif- 
ferent valence occurring upon full CO. (c) The polarization 
Px has the same size as the a Bravais vector with 1 electron 
charge (27r in Berry phase), so that Px is zero; on the other 
hand, the net Py is non-zero. 
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FIG. 4: Perspective view of 00 of ¥e-t2g down-spin states 
(within an energy range up to 1 eV below Pf) at (a) COII 
and (b) COIII. The strong/weak next-neighbor orbital over- 
lap and the consequent intersite Coulomb repulsion is high- 
lighted. 



t2g-t2g ovcrlap. Figure [4] shows the three dimensional net- 
work of Fe-t2g 00. It is clear that the t2^-orbital relative 
to the "extra charge" on top of Fe^+ site at COII pattern 
is aligned along a diagonal direction in the octahedron, 
so as to avoid the strong inter-site Coulomb interaction 
with neighboring t2^-orbitals at Fe^+ sites. Therefore, we 
argue that the COII pattern is energetically lower than 
COIII, due to minimization of Coulomb repulsion. 

CO-induced Ferroelasticity - From the above dis- 
cussions, we conclude that the ground state of 
Ko.6Fe[)^gFe[)^iF3 is FE with a polar COII pattern (al- 
though this has yet to be confirmed experimentally). Fur- 
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thermore, the polar CO patterns induce "unbalanced" 
(from the valence point of view) planes at z = 1/4 and 
at z = 3/4, so that a monoclinic crystal distortion is ex- 
pected (in this case along the direction of P, which is de- 
termined by the COII pattern as well). This is proposed 
as responsible for the reported ferroelastic phase Q. To 
check this effect, we optimized the a, 6, c lattice vectors 
under CO-I and CO-II patterns, keeping the cell vol- 
ume fixed to the experimental value. It turns out that 
the optimized lattice is in fact monoclinically distorted, 
with an angle Z6c=90-0.059° at CO-II, and an Zac=90- 
0.040° at CO-III. This reveals a significant fact: ferroe- 
lasticity is strongly coupled with CO, which in turn de- 
termines the direction of the FE polarization. We note 
that the ferroelastic distortion in the ac or he plane is 
peculiar to this CO system and distinct from the con- 
ventional monoclinic distortion in the ah plane, the lat- 
ter occurring in other ferroelastic TTB compounds, e.g. 
Ba2NaNb50i5, as well. Although these monoclinic dis- 
tortions in Ko.6FeQ^gFeQ^4F3were not experimentally ob- 
served, we remark that they are rather small and possibly 
below the detection limit. 

Ferroelectric anisotropy - The proposed polar COII 
pattern may not be spontaneously long-ranged; how- 
ever, it may be possible to stabilize large FE domains by 
field-cooling upon applying Ey, in analogy with LuFe2 04 
(where the FE and the anti-FE states are energetically 
close and the stabilization of the FE phase occurs via 
field-cooling). We further speculate that the energetically 
competing CO-III pattern inducing a large Px might be 
realized by applying E^. In the case of a (strong enough) 
applied electric field rotating in the xy plane, the induced 
P along the field would therefore change its saturation 
value. This "ferroelectric anisotropy" (FEA), denoted 
as the energy required to modify the direction (as well 



as size) of the permanent polarization by switching the 
crystal between different CO phases (in this case COII 
and com), may found applications in future devices, 
such as multiple-state memories where the information 
can be stored by exploiting not only the sign of P, but 
also its direction. The FEA is peculiar for "improper" 
ferroelectricity induced by electronic degrees of freedom; 
in analogy with CO-induced polarization, the possibility 
to control the direction of P in spin-spirals manganites by 
means of a magnetic field was already proven and sugges- 
tions towards devices harnessing FEA already came. [2^ 

In summary, tungsten-bronze systems, previously 
known for hybridization-driven ferroelectricity occur- 
ring in Ba2NaNb5 0i5, branch into the class of improper 
multiferroics. Indeed, we have put forward CO as the 
origin of ferroelectricity (as well as of some ferroelastic 
modes) in a non-oxide TTB compound, Ko.6FeQ^gFeQ^4F3. 
Several energetically-competing CO patterns were pre- 
dicted from first-principles, with potentially different 
directions and magnitudes of the polarization. This 
"ferroelectric anisotropy", typical of electronic mag- 
netic ferroelectrics, shows a high technological appeal 
in multiple-state devices. In addition, the strong inter- 
play between CO, ferroelasticity and ferroelectricity in 
Ko.6Feo^gFeQ^4F3 makes it an excellent compound where 
multiferroic effects are manifestedly at play. 
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